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Abstract 

We propose a new multi-component two-dimensional Toda lattice hierarchy (mc2dTLH) which in- 
cludes two-dimensional Toda lattice equation with self-consistent sources (2dTLSCS) as the first non- 
trivial equation. The Lax representations for this mc2dTLH are presented. We construct a non-auto- 
Backlund Darboux transformation (DT) for 2dTLSCS by applying the method of variation of constant 
(MVC) to ordinary DT of 2dTLSCS. This non-auto-Backlund DT enables us to obtain various solutions 
such as solitons, rational solutions etc., to 2dTLSCS. 

1 Introduction 

Multi-component generalizations of soliton equations attract a lot of attention from both physical and math- 
ematical points of view [1-8]. The multi-component KP (mcKP) hierarchy given by [1] contains physically 
relevent nonlinear integrale systems such as Davey-Stewartson equation, two-dimensional Toda lattice and 
three-wave resonant interaction equation. The multi-component Toda lattice Hierarchy [8] contains non- 
abelian Toda lattice equation. There exist several equivalent formulations of this multi-component soliton 
equations. For example, there are matrix pseudo-differential operator (Sato) formulation, r-function ap- 
proach via matrix Hirota bilinear identities, multi-component free fermion formulation for mcKP hierarchy. 
For two dimensional Toda lattice hierarchy (2dTLH), a similar matrix-difference operator approach to multi- 
component hierarchy was also presented by [8] . 

Another kind of multi-component generalizations to soliton equations are the so-called soliton equation 
with self-consistent sources (SESCS), which were initiated by V.K. Mel'nikov [9-11]. For two dimensional 
Toda lattice equation (2dTL), the corresponding 2dTL with self-consistent sources (2dTLSCS) was first 
presented in [12] by source generating method as follows 
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In [13], we proposed a method to construct a new multi-component KP hierarchy which includes two 
kinds of KP equation with self-consistent sources presented by Mel'nikov [9]. In this paper, we will present 
a new multi-component 2dTLII which includes 2dTLSCS as the first non-trivial equation. 

We briefly recall the framework of [8] as follows. Let x = (xi, 2:2, • ■ • ) and y = (yi, 2/2, • ■ • ) be two series 
of variables. Let n G Z be a discrete variable. Then Lax equation of 2dTLII is given by 

L,„, = [B^,Ll (2a) 
Ly^^[C^,L], (2b) 
M,„ = [B„,M], (2c) 
My^^[C^,M], (2d) 

where 

L = A + UO + uiA^^ + U2A"^ H , 

M = v^iA^^ +V0 + viA + V2A^ H , 

A is a shift operator such that Af{n) — f{n + 1)A, Ui and Vi are functions of x, y and n, B,n — L™ stands 
for the positive part (> 0) of with respect to the powers of A and Cm = stands for the negative 
part (< 0) of M™. The commutativety of ^ gives rise to zero-curvature equations of 2dTLII 

Bk,x„^ — Bm.Xk + [Bk,Bm] = 0, (3a) 
~ Cm.yk + [Cfe, C,„] — 0, (3b) 
Bk,y^ — Cm,Xk + [Bk,Cm] = 0. (3c) 

When m = k ^ 1, ^ leads to following 2-dimensional Toda equation 

Uy = V — V , (4a) 

/ (-1)- 

Vx ^ V {U - I 

and Lax pair for Q is 



Vx V [u - u ) , (4b) 



^P, = B{^) = iA + u){^), (5a) 
ijy = C{^) = {vA-'m, (5b) 

where x :— xi, y :— j/i, B :— Bi, C :— Ci, u :— uq, v :— Eliminating u from ([4]) and introducing 

q := q(n, x, y) which satisfies 

V := exp {q - q , (6) 
then ^ gives the so-called two dimensional Toda lattice equation: 



exp 



{q-q - exp {q^^^ - q^ . (7) 



Our multi-component generalization to 2dTLII can be presented as follows. We first introduce a new 
vector field dy^ , which is a linear combination of all vector fields dy^ . Then we get a new Lax type equation 
which consists of dy^-^iow and evolutions of wave functions. Under the evolutions of wave functions, the 
commutativety of dy,^, dy^ and d^^. flow give rise to the new multi-component 2dTL hierarchy (mc2dTLH). 
This hierarchy enables us to derive the 2dTLSCS in different way from [12, 14] and to obtain their Lax 
representations. This hierarchy is also different from mc2dTLH given by [8]. 

In the second part of our paper, we solve 2dTLSCS by means of Darboux transformations (DT). Since 
the Lax representation for 2dTLSCS is obtained, we can construct an auto-Backlund DT for 2dTLSCS, 
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which transforms between the solution of 2dTLSCS with same number of source terms. However, such auto- 
Backlund transformation can not be used to construct non-trivial solution from the trivial solution. The 
idea for us to construct non-auto-Backlund DTs is to consider 2dTLSCS as 2dTL with non-homogeneous 
terms (i.e. self-consistent source terms). Inspired by ODE method, we can apply the method of variation of 
constant (MVC) to DT to find a new non-auto-Backlund DT which transforms original solution of 2dTLSCS 
with N self-consistent sources to a new solution of 2dTLSCS with {N + 1) self-consistent sources. This 
non-auto-Backlund DT enables us to find various solutions for 2dTLSCS. Furthermore, we obtained the 
m-time repeated non-auto-Backlund DTs formula, and exhibit some solutions of 2dTLSCS which include 
solitons, rational solutions and etc. 

Our paper will organized as follows. In section [2] we propose resolvent identities and present new 
mc2dTLH which includes 2dTLSCS. In section [3] we first construct auto-Backlund DT for 2dTLSCS. Then 
by applying method of variation of constant to this DT, we find a non-auto-Backlund DT which can increase 
the number of source term by 1. We obtain m-time repeated non-auto-Backlund DT formula which can 
be expressed in compact form with Casoratian determinants. In section |4l we present some solutions to 
2dTLSCS by using this m-time repeated DTs. 

2 New Multi-component 2-dimensional Toda Lattice Hierarchy 

2.1 Sato approach and resolvent identities 

First we introduce some useful notations and definitions which can be found in [8]. 
Definition 1 (The residue of shift operator). Let P = X]iez^«^*' then residue of P is 

RcsaP Po- 

Definition 2 (The adjoint operator *). P* = J^tei.^^'^^- 

Definition 3 (Shift operator's action). The shift operator action P(A") can be defined as 

Definition 4 (Formal inversion of difference operator A). For difference operator A = A — 1 , the formal 
inversion are given by 

A^^^-J2^' or AI^ = ^ A\ 

i>0 i<-i 

Introduce wave operators 

^ bo + biA-^ + b2A-^ + ■ ■ ■ , 
=co + ciA -I- caA^ H , 

where bo = 1. Ueno and Takasaki proved [8] that if L and M are solutions to ^ then there exist wave 
operators, such that L and M can be written as 

L = W AW \ M = W A'^W (8a) 

and 

dy^w'^' = M^oW'^', dy^w"" = -M^oW"" ■ (8c) 
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Define wave function 

j:>0 



i>0 

and adjoint wave function 

where ^(x, A) — X]i>i s^iA', ^(y, A^^) = X]i>i Vi^^^ ■ Then ([2]) can also be given as the compatibihty condition 
of the following linear evolution equations 

Lw — \w , Mw — \ w , (9a) 

dj;^^ = BmW , Ca;„it; ^ -B,nW , (9bj 

Lemma 1 (Ueno and Takasaki [8]). Suppose P — ^P^A*, Q — '^Qjh!' , then 

ResA P Q* = RcsA \-^P{X') ■ Q(A-"). 

Proof. Only need to show for P = PiA\ Q = Q^A^. ResA PQ* ^ Rcsa PjA'-^Qj = S,,jPiQj, while 
ResA A-ip(A")Q(A-") = Res^ A-iF,A"+*gjA-"-^' = 6,^jP,Qj. □ 

Similar to the KP theory, in which the principle part of resolvent can be expressed in terms of a quadratic 
form of wave function and adjoint wave function [15], we have the following resolvent identities for 2dTLH. 

Proposition 1 (Resolvent identities). 

k>o fcez 

keZ k>0 

Proof. We prove one of them, others are similar. Since L — W' 'AW'' 

m>l 



J2 ResA [w'^'A'^w'^'-'A"^'j A-™ 

rn>l 

J2 Rosa A-1 A'=(A")e«(-^^)) • (a-™!^*""' *-i(A-")e-«-'^)) A" 



m>l 

EResAA'=-iw;'°°'A-™w'°''' 

m>l 

ResA A'=-1w*°°'ai1u.'°°'*. 
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2.2 New mc2dTL hierarchy 

For fixed k > 1, N > 0, we define a new time variable yk sucli that the corresponding vector field is 

N 

dy.^dy,+Y,Y.^9y^ (10) 

1=1 j>l 

where Xi are distinct arbitrary non-zero parameters. Then the ijk flow is given by 

—L^[Ck,L], M^[Ck,M], (11) 

with 

N 

i=l j>l 

which, according to Proposition 1, can be rewritten as 

N 



(0) ^ (0)* 



By setting Wi = , w* = , the compatibilities of Q and (jlip give rise to the following new multi- 
component two dimensional Toda lattice hierarchy. 



Proposition 2. We have the following new mc2dTLH, 
for m ^ k: 



for m = k: 





— Bk,x„, + [Bm, Bk] 


= 0, 








c 


~ (^k,ym + [C'mi Ck] - 


= 0, 










— Ck,xm + [Bm, Ck] 


= 0, 










~ Cm,Xk + [Bk, Cm] 


= 0, 










^ Bm[Wi) Wi,yr„ 


— Cm( 


'Wi) {i = 


1,.. 


■,N), 




= -Bmiw*), W* 


,ym ^ " 


-CmiO, 








- Ck,Xk + [Bk,Ck] — 













= Bk{w{), dx,w* 


= -Bl 




= 1,. 


..,N) 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 
(12f) 



(13a) 
(13b) 



N 



where Cfe = + ^ WiA_^ 



It is worth noting that Wi and w* need not necessarily to be the wave function and adjoint wave function. 
In fact, the equations (|12ep and (I12f|) ( or (|13bp ') ensure the closeness of (I12ap - (ll2d[) (or (|13ap ). Furthermore, 
under the conditions (|12ep and (|12f|) (or (jl3bp ). one can easily obtains the Lax representations of (|12ap - (jl2dp 
as 

ijx^=Bm{ip), iJx,=Bki^p), (14a) 

i^Vrr.^ Cm ^y,=CkW, (Ub) 

or get the Lax representation of (|13ap as 

^x,=Bk{i^), ^Py.^Ckii^). (15) 
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Example 1 (Two dimensional Toda lattice equation with Self-consistent sources). When m = k = 1 , let 
u = uo, V = v-i, X = xi, y = yi 

Bi^A + u, Ci = vA-\ 

then becomes 



N 

Uy = —A{v + ''^WiW* ), Vx = v{u - u ), (16a) 



w,^,^B,iw,), wl^ = -Bl{w*), i = l,...,N (16b) 



Under u — Qx, v — exp(g — '), U6]) yields 

N 

qxy = e"-"''"' -e9'"-9 + ^(u;,u;*),, (17a) 

i=l 

Wi^x ^ + QxWi, (i = l,...,iV) (17b) 
wlx = " - qxw*. (17c) 

This is two dimensional Toda lattice equation with N self- consistent sources (2dTLSCS). 
Analogously, let us introduce x^, such that 

N 

dxk = dx^ + ^^^J^dx,, 

i=i j>i 

then we will get another new multi-component two dimensional Toda lattice hierarchy. 

Proposition 3. We have another new mc2dTLH as follows, 
for m ^ k 

(18a) 
(18b) 
(18c) 
(18d) 

N, (18e) 

(18f) 



for m = k 



where 



^rn.Xj^ 


— Bk,x„, 


+ [Bm, Bk] 


= 0, 




Cm,Xk 




+ [Cm, Bk] 


= 0, 






~ ^k.ym 


+ [Cm, Ck] 


= 0, 








+ [B,n, Ck] 


= 0, 






= Cm(m, 




= Bm 


(wi), i 




= -C™C 






-Bm{w*) 


Ck,Xk 




-[Ck,Bk] = 


0, 





(19a) 

dy,w, = Ck{wi), dy,w* = -Cl{w*). i = l,...,N, (19b) 



N 

Bk^ Bk-^w,l\^^w*. 

i=l 

Example 2 (2dTLSCS [12, 14]). When m = k = I, (T^ leads to (Qp. It is interesting to see that is 
equivalent to ( [_?7| j under 

x-> -y, y ^ -X, q~> q, 
Wi —e'^w*, w* — > e^'^Wi. 

So hereafter we may concentrate on 2dTLSCS U6]) . This transformation was discovered by Prof. Hu Xing- 
biao. 
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3 Darboux transformation for 2dTLSCS 



In the second part of our paper, we concentrate on 2dTLSCS (|16p. First recall the Lax pair of 2dTL equation 
([5]). For convenience, hereafter we denote B — Bi, C = Ci, u = uq, v — u-i, x = xi, y = yi. 

3.1 Applying the method of variation of constant to DT of 2dTLSCS 

Let us first introduce the notion of Casoratian determinant: for m discrete variable functions hi, ■ ■ ■ ,hm, 
the Casoratian determinant 

hi ■ ■ ■ hjn 

(1) , (1) 



cas(/ii, • ■ • 



h\ 



Darboux transformation for 2dTL (j4|) was given in [16]. Let us first recall this DT and its proof as 
following Lemma. 



Lemma 2. Let h be special solution to (0). Let P = A + a.cr := —h^ ^ /h, then DT 

u := u + (7 — a , (2(Jaj 

i:^va/a''\ (20b) 

7 >^ cas(/i, ip) , 

^ := V{ij) = -^j^, (20c) 

gives a new solution to Thus u, v are new solution for Q). 

Proof. Since B := A + u, C := vA^^, t/j = T>{t/j), a sufficient condition such that ^ holds is 

n^+VB- BV = 0, (21a) 
Vy+VC -CV = 0. (21b) 

Notice that 

V{h) = 0, (22) 

take partial derivative dx,dy to ([22]) . one gets 

V,{h) + V{h,) = (V, + VB){h) = BV{h) = 0, 

Vy{h) + V{hy) = [Vy + PC)(/l) ^ CV{h) = 0, 

which mean 

{Vx +VB - BV) (h) ^0, (23a) 
{Vy + PC - CV) (h) = 0. (23b) 

From (|20ap and (j20bp one knows the operators acting on h in are scalar functions multiplications. So 
(I2TI1 holds. □ 



The Lax representation for 2dTLSCS ([16]) is 

V'x = S(V), (24a) 

i^y = {C + J2^Az'w*m. (24b) 
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Note that Lax representation holds under following equations 

w.,^,=B{w,), i = l,...,iV, (24c) 
wl^ = -B*{w,). (24d) 

Proposition 4 (Darboux transformation for 2dTLSCS (|16[)). Let h be a special solution to {24% V = \ + a , 
a := —ft,' '//i. Based on the Darboux transformation 12(J\). define 

Wi := V{wi) = , (25a) 

h 

^l:^V*-\K)^-^-^, (25b) 

where S :— AAl^. Then i20\) and h25]) together give a new solution to ^24\ )- Thus one gets a new solution 
to (01). 

Proof. From Lemma [2] it is easy to see that Wi defined by p5ap satisfies (|24c|) . It is necessary to prove w* 
satisfies (I24dp . From the proof of Lemma [2] we know 

(9, - B)V = V{d, - B). 

Taking formal adjoint * to this equality and rewrite it as 

{-d^ - B*)V*^^ ^ V*~^{-d^ - B*). 

This is a sufficient condition for T>*^^ to be the Darboux transformation for (j24dp . Thus we have proved 
(j25bll . At last we need to prove that Darboux transformation given by ((20|) and ([25| fulfills (j24bp . That is 

N N 

Vy+VC + J2 'DmAZ^w* -CV~Y^ mAZ^w*V = 0. (26) 

i=l i=l 

Based on (j21ap . we have to prove the extra terms w.r.t. Wi, w* in (|26p are equal. For every i, we have 

(1) Q/T *X , (1) 

- ^ + -^u;*A_i(/iO + AAjw* 

—WiAZ^w* - w,AZ^w*A + w,AZ^w*—r- 

h h 

, (1) 

' A -1 * ~ A -1 S(hw*) , . . _i S(hw*) 

h h h 

, (1) 

= - ^ S(ft<) + u;f ' AA:i< - —w.KAZ^w* 

+ w,A_K ^ ~ WiA_ \ - w^AAjw* = 0. 
h h 

□ 

Theorem 1 (Darboux transformation and method of variation of constant for 2dTLSCS ([U])). Let f and g 

he two linear independent solutions to |^^[ ). Suppose a(y) is arbitrary functions of time y. Let h f + o-iy)g, 

WN+i = c{y)V{f), (27a) 

d{y) s 

^N+i = (27b) 
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then i2U\) . \25]] and \27\j give a new solution for H6\) and \24\j with TV + 1 self- consistent sources , where 
c{y), d{y) satisfy c{y)d{y) = dy\oga{y). 

Proof. It is easy to see wn+i satisfies (I24cp . To prove that Wtv+i satisfies (|24d[) . we have 



^(1)2 



(h +u h ), 



d{y) (1) _ h h_ d{y) _ d{y)h _ (d d{y) 



Based on proposition 31 we want to show that extra terms come out from 



V 



can be cancled out. It because 



(1) , (1) 

h 



WN+iA_^iv*j^_^_i(A - h /h) 



=^ cas(5, M - c(2/My)^^^ AI^/."' -^A - i) 
=^ cas(5, /) - c{y)d{y)a{yf-^^^ = 0. 



□ 



3.2 m-time repeated non-auto-Backlund DTs 

Theorem 2. Let fj and gj (j ^ 1,2, . . . ,m) be m pairs of independent solutions to {^p. Suppose aj{y) are 
arbitrary functions of time. Let 

h.j := fj+aj{y)gj. 

Then after m-time repetition of Theorem]^ we find a solution for {24^ with N -\- m self- consistent sources, 
which is 



where 



cas(/ii 



v[m 
Wi [m 
w* [m 
WN+j[m 



1 ) 



(m) cas {hi,---,hjn) cas{hi,- ■ ■ ,hm) 



cas'^' [hi 

(1)/, , N (-1), 



tl, • ■ ■ , ^m) Cas(/li, • • • , hm) 

cas'"' (^1, ■ ■ • , hm) cas' {hi, ■ ■ ■ , km) 



cas2(/ii, • • • , hm) 
cas(/ii, • • ■ ,hm,Wi) 



i = 1 



cas(ft.i, • • • ,hm) 

cas(/ii, • • ■ , hm, w*) 



, . . . , ^ . , 



N, 



(-1) 



cas'"(/ii,-- - ,hm) 



cas(/ii, • • ■ ,h„ 



{~l)^-^d,{y) 



hi 



(™-2) 



cas [hi. 



,hj,- 



) hm) 



hi 

hT 



cas' '{hi-- - ,hm) 



,cas(/ii, • -- ,hm, f) 



S{hif) 
hi 



S{hmf) 

, (1) 



(28a) 

(28b) 
(28c) 
(28d) 
(28e) 
(28f) 
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and Cj{y)dj{y) = dy\ogaj{y). 

Proof. Since each time Darboux transformation has the form T) = A + a, after m-time repetition, corre- 
sponding operator has the form 



2?(m) = A'^ + a^^iA" 



+ Co- 



There are m indetermined coefficients ai, i — 0, . . . ,m — 1. From (|22p we know 



So the indertermined coefficients satisfies 



hi h'l ' 

(1) 



h2 /l2 



By Cramer rule, it is easy to see 



i = 


= 1,2,. 


. , m. 










hi 




CTl 












, 



ao = (-1)' 



, cas (hi, - ■■ ,hm) 
cas(/ii, • • • , 



cas(fei, ■ ■ ■ , hjn) 
cas(/ii, • • • , h,n) ' 



Note that DT of u, v are the same for 2dTL and 2dTLSCS. So we may omit source term temporary. That 
is, assuming that P(m) transforms ^ to — ^[m], and satisfying ■i/'j. = B^p and — Ctp, then we have 



V{m)^ + V{m)B - BV{m) = 0, 
V{m)y + V{m)C - C'P(m) = 0. 



(29a) 
(29b) 



Comparing the coefficient of A™ in (|29ap . we have (|28ap . Comparing the coefficient of A^^ in (j29bp . we 
have (j28bp . For arbitrary eigenfunction w, its DT w — 'D{m){w) can be expressed in a compact form 
(j28cp according to the Laplace expansion formula. For (j28dp . we need induction. Suppose for any adjoint 
eigenfunction iti*, the m-time DT formula is correct, then by (|25bp . the m + 1-th DT is 



w* [m + 1] 
S 



S{h^+i[m]w*[rn]) 



/im+i[m] 



(1) 



hm+i[m 



(1) 



hm+i[m- if^ S(/i„i[to- l]w*[m- 1]) S(ft.m+i[m- l]w*[m~ 1]) 



hm+i[m 



(1) 



By assumption 



S{hm[m - l]w*[m - 1]) _ ^ j^-^m+i ca5(/ti, ■ ■ ■ ,hm,w*) 
/im[r7i-l]"' cas'''(/ii, • • • ,/im) ' 

S(ftm-n['-» - l]w*[m - 1]) _ ^ ^^,„^i Cas(/ll,--- ,/i-m-l,/tm+l,W*) 

/im+i[TO-l]"' cas"'(fti,--- ,/i,„_i,/i„i+i) 



we know 



w*[m + 1] 



(-1) 



m+l 



CaS (/ll, • • ■ , /im-l) CaS {hi, ■ • ■ , /im+l) 
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cas(/ii, • ■ • , ft-m, w*) cas {hi,- 



, ilm-l, "■m+l J 



- cas(ft.i, • • ■ , hm-1, hm+i,w*) cas {hi, • ■ • , /i^) 
cas*'' (/ii, • • • , cas*'' {hi, ■■ ■ , h^+i) 



S{hiw*) 
hi 



hr 



S{hm~lW*) S{hmW*) S{hm+lW*) 



-1) 
-1 



(™-l) 



(1) 



hi 



-1) 



, (™) 



'm+1 



(m-l) 

^m+1 
, (1) 
"m+1 





(1) 



m+l Cas(/j.i, 

CaS*''(/li, • • • , ft-m+l) 



Next we want to prove the m-th repetition formula for new source terms. Firstly, it is easy to see 
WN+j[iTi] = Cj{y)fj[m] can be derived from (|28ep . For w'^_^_j[m], we use induction. Suppose when j < m the 
formula for w^^^ [m] is given by (j28f|l . then 



JN+j 



[to + 1] = 



/lm+l[TO]*'' 



Because hm+i[m\ is obtained by m-time repetition of DT by using hi ■ ■ ■ , h„i sequentially, which is equivalent 
to TO-time repetition of DT by successively using hi, ■ ■ ■ , ft-j-i, ft-j+i, • • • , h^ and at last hj, 

u r 1 J, r ii<" hj[m~lf'\ 

hm+l [m\ = hm+l [to - IJ —hra+l [to - IJ . 



hj [to — 1] 



Note that (|28fj) . we have w 



N+jy 



hj[m-lY 



wli^j [to + 1] = - 



do{v) 



hm+i[m] 



SA 



hj [to - 1] 



(- -^^m+l-.7- '^3(?/)caS {hi,--- ,hj,--- ,hm+l) 
cas*''(/ii, • • • , hm+l) 



When j = TO + 1 , 



+m+ 



i[m+l] 



dm+i{y) d,n+i{y) ca.s {hi,---,hm) 



hm+i[m] 



cas*''(/ii,--- 



□ 



4 Solutions for 2dTLSCS 

Let us start from trivial solution (7=l,u = l,u = 0, A'' = Ofor 2dTLSCS (|24p . The Lax pair reads 



(-1) 



(30a) 
(30b) 
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4.1 Solitons 

Equations (j30p have two linearly independent solutions 

/(n, X, y) — exp{nuj + zx + z^^y), g{n, x, y) — exp(--?iijj + z^^x + zj/), 
where z = e". Let a(y) = e"^^\ then 

h = f + a(y)g = 2 exp • cosh Z, 

here 

f2 = coshw • X + cosh a; • y + a/2, Z = nuj + sinhw • x — sinhcj • y — a/2. 
By (pS)) . taking m = 1 we have the following 1-soliton solution for 

cosh(Z + 2w) cosh(Z + 



cosh(Z + uj) cosh Z 

cosh(Z + uj) cosh(Z — uj) 



w*[l] 



cosh^ Z 
sinh • 
cosh Z ' 



2cosh(Z + cj)' 



where c{y)d{y) — a. 

If take two pairs of independent solutions, with respect to Zj = e'^^ {j = 1,2), i.e. 

fj = exp{nujj + ZjX + zj^y), gj = exp(-na)j + zj^x + zjy) j = 1, 2. 



Let Qjiy) = e°'^'-y\ then 
where 



hj = fj + ajgj = 2exprij • coshZj, 



flj — coshwj • X + coshwj • y + Q!j/2, Zj ~ nujj + sinhojj • x — sinh a; j ■ y — aj 
To simplify the notion, for fc G Z, define 



Hk = 



cosh Zi cosh Z2 

cosh(Zi + kuti) cosh(Z2 + kuj2) 



= sinh 



k{uJi - UJ2) 



sinh I Zi + Z2 + -{uji + LJ2) 



sinh 



k{uji + LJ2) . 



sinh [ Zi - Z2 + -{uji - UJ2, 
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Then 2-soliton solution for (fTC)) is 

u[2] 

v[2] 





H2 




Hi 




-1) 





, ^ , , 2sinha;i(coshwi — coshw2) expili cosh(Z2 + 
u>i[2] = ci{y)ai{y) 



W2[^] = C2(y)a2{y) 
wl[2]-- 



2 sinh UJ2 (cosh oji — cosh L02 ) exp $12 cosh(Zi + uji ) 



Hi 

di{y) exp(-f2i) cosh(Z2 + CJ2) 



^2 [2] = 
where Cj{y)dj{y) — aj. 



^2(2/) exp(-il2) cosh(Zi +UJ1) 



4.2 Rational solution 

In equation (j30p . noticing that 9^?/; is another solution. Since x, y) = z" exp(za;+z^^j/) and /fc(n, a;, y) := 
d^g {k > 1) are all independent solutions for ([501) . Let ^ := zx + z^^y, then 

fi{n,x,y) = = z""^e^(n + z^^), 

/2(n, a;, j/) = d'^g = z"-2gC(„2 _^ ^nz^, -n + z^il + z^^^^), 
f3{n,x,y) = • • • 

Let h^. — ft; + a{y)g. Take /c = 1, m = 1 in (j^S]) . one yields 

z 

"^^^ " "(7; + za+l/2)2-l/4' 
v[l] = l- 



(?7 + za) 



z"+ie«a 



7; + za 



77 + za + 1 ' 

where r] = n + z^z,ciy)d{y) — logo. This is a rational solution for ([TB)) . 
If take fc = 2, m = 1, we find another rational solution for 



r,, V + z 0, ri + z^a 

u[l\ :— z ' 



rf ' + z^a 1] + z-^a 

(t] ' + z^a^ij] ' + z^a) 
(ji + z^a)^ ' 



w[l] 2c(y)az«+ie«^^±4^, 
77 + z^a 

7] + z^a 

where 77 = ri^ + 2nz^2 -n + z^^^ + z^^^^, c{y)d{y) = ^ log a. 
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4.3 Other solutions 



Let 

J — ^n^zx+z^^y ._ ^n^F{x,y,z) ^ _ ^-n^z^^x+zy ._ ^-n^G(x,y,z) 

be pair of solutions to (|5D1) . then and are another pair of solutions to ([501) . Let 

h — f + a{y)g = 2 exp CI ■ cosh Z, 
where and Z are defined in subsection 14. II Then 

hz = fz+ a{y)gz = 2f7^e" cosh Z + 2e^^Z^ sinh Z. 
From (|28p . taking m — 2, we construct solutions with singularities. For simplicity, we define 



h hz 

,(fc) 



fc k 

ie'^^Ck, whereCfc = {Zz-\ ) sinh(fcw) H sinh(2Z + kuj). 

2z 2z 



where 



sinh^ u! 

cas(/i, hz, f) = -8a{x)e^^^ cosh(Z + uj), 

z 

cas{h, hzjz) = 4a(a;)e^" (^d['^ coshZ + Di cosh(Z + 2oj) - D2 cosh(Z 

u2 



,n k ^ ,n 

Dk = -(- H h H Gz) sinh(fcw) H 5- sinh(A;w) H cosh(fcw), fc = 1, 2. 

z 2z z 2z Az z 

then find the following solution for 



„(i) „ ^(i)^<-i) 
w;i[2] = -2ci{y)a{y) — cosh(Z + w), 



( Dl coshZ + Di cosh(Z + 2w) - D2 cosh(Z + cj) 

W2[2] = C2{y)a{y)- 



wl[2] = -d,{y) 
wl[2]^d2{y) 



Ci 

cosh(Z + tj) + (Zz + 1/z) sinh(Z + w) 



2e^&^' 

cosh(Z + f^) 



where Cj{y)dj{y) = ^loga(y). 

Conclusion 

We present a new multi-component two dimensional Toda lattice hierarchy, which enables us to find the 
two dimensional Toda lattice equation with self-consistent sources in different way from [9-12,14] as well as 
their Lax representations. Since two dimensional Toda lattice equation with self-consistent sources can be 
considered as a two dimensional Toda lattice equation with non-homogeneous terms, method of variation 
of constant can be applied to the ordinary Darboux transformations for 2dTLSCS to construct a non-auto- 
Backlund Darboux transformations. Then it offers a different way to solve 2dTLSCS in contrast with [12,14]. 

The 2dTLH offers various types of reductions, for example periodic reductions and reductions to Toda 
lattice equation. It is an interesting question does this new mc2dTLH offers similar reductions. We may 
discuss such problems elsewhere. 
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